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A  REDUCED-ORDER  MODEL  OF  A  FLUID-LOADED, 
ELASTOMERIC-COATED  METALLIC  PLATE 


1.  INTRODUCTION 


Elastomeric-coated  metallic  plates  are  used  in  a  wide  variety  of  applications.  In  marine 
systems,  they  are  typically  used  to  absorb  energy  or  to  transmit  acoustic  energy  from  the  exterior 
of  an  underwater  object  to  a  location  on  the  interior  of  the  coating  or  metallic  plate.  Sometimes 
sensors  are  embedded  in  the  elastomer  or  the  metallic  plate  to  measure  the  energy  field  that  is 
insonifying  the  stmcture.  Energy  at  different  frequencies  and  wavenumbers  will  cause  a 
different  response  in  both  the  coating  and  the  metal.  One  way  to  model  this  response  is  a  full- 
order  model — i.e.,  deriving  the  equation  of  motion  of  two  thick  elastic  plates  and  then  coupling 
them  to  an  acoustic  wave  equation  that  models  the  fluid  load.  Because  the  resulting  system 
matrix  is  eight  by  eight,  this  model  has  a  final  solution  that  is  expressed  in  open  form,  and  the 
system  matrix  must  be  numerically  inverted  at  every  frequency  and  wavenumber  of  interest.  It  is 
difficult  to  understand  the  physics  in  such  a  solution  because  the  interaction  of  the  coefficients  in 
the  response  is  implicit.  A  second  method  is  to  model  the  metallic  plate  as  a  discrete  mass, 
which  eliminates  one-half  of  the  equations.  Once  this  is  done,  the  system  equations  have  a 
dimension  of  four  by  four,  which  allows  the  matrix  to  be  symbolically  inverted,  yielding  a 
closed-form  solution  to  the  problem  at  every  wavenumber  and  frequency.  Closed-form  solutions 
are  advantageous  because  the  dynamics  associated  with  the  system  are  explicit.  In  addition, 
applying  inverse  methods  to  the  problem  is  much  easier  when  a  closed- form  model  is  available. 

This  report  develops  both  a  full-order  model  and  a  reduced-order  model  of  a  fluid-loaded, 
elastomeric-coated  metallic  plate  and  compares  the  results  graphically  for  system  displacements 
and  stresses  at  various  frequencies  and  wavenumbers. 
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2.  FULL-ORDER  MODEL 


The  transfer  functions  are  derived  by  modeling  the  elastomeric  coating  and  the  metallic  plate 
as  elastic  structures,  both  of  which  are  independently  governed  by  the  equation* 

//V**u+(;i+//)vv«u  =  /7^,  (1) 

a 


where  p  is  the  density,  X  and  p  are  the  complex  Lame  constants,  t  is  time,  •  denotes  a  vector  dot 
product,  and  u  is  the  Cartesian  coordinate  displacement  vector  expressed  as 


u  = 


'u^{x,y,z,t) 
Uy{x,y,z,t)  L 
u.{x,y,z,t) 


(2) 


with  subscript  x  denoting  the  direction  parallel  to  the  plates,  y  denoting  the  direction  orthogonal 
to  the  x-direction,  and  z  denoting  the  direction  normal  to  the  plates.  The  modeled  geometry  and 
the  coordinate  system  of  the  structure  are  shown  in  figure  1 .  The  symbol  V  is  the  gradient  vector 
differential  operator  written  in  three-dimensional  Cartesian  coordinates  as^ 


d  .  d  .  d  . 
V  —  —  l  H - /  H - Z, 


(3) 


with  ix  denoting  the  unit  vector  in  the  x-direction,  iy  denoting  the  unit  vector  in  they-direction, 

and  4  denoting  the  unit  vector  in  the  z-direction.  The  symbol  is  the  three-dimensional 
Laplace  operator  operating  on  vector  u  as 

V^u  =  +^^Uyiy  +  ,  (4) 


with  operating  on  scalar  u  as 


.  d^u 


x,y\z 


dyr 


dz^ 


(5) 
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The  term  V»u  is  the  divergence,  which  is  equal  to 


V«u  = 


^  dy 


+ 


du. 

dz 


Figure  1.  Modeled  Geometry  and  Coordinate  System  for  Full-Order  Model 


(6) 
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The  applied  loading  effect  of  the  pressure  in  the  fluid  acting  on  the  surface  of  the  top 
elastomeric  plate  is  modeled  as  an  acoustic  load.  Note  that  p.  A,  and  p  are  the  properties  of  each 
specific  layer. 

The  displacement  vector  u  is  written  as 

u  =  V^  +  Vx^?,  (7) 

where  is  a  dilatational  scalar  potential,  x  denotes  a  vector  cross  product,  and  y/  is  an 
equivoluminal  vector  potential  expressed  as 


y/  = 


y/y{x,y,z,t) 

[y/.{x,y,z,t)] 


(8) 


The  problem  is  now  simplified  as  a  two-dimensional  response  (j  =  0  and  di^^l  dy  =  0), 
since  it  is  known  that  an  acoustic  plane  wave  excites  the  plate.  Expanding  equation  (7)  and 
breaking  the  displacement  vector  into  its  individual  nonzero  terms  yields 


w^(x,z,0 


d^{x,z,t)  dy/y{x,z,t) 
dx  dz 


(9) 


and 


dz  dx 


(10) 


Equations  (9)  and  (10)  are  next  inserted  into  equation  (1),  which  results  in 


c]^^(f){x,z,t) 


d^(f){x,z,t) 


(11) 


and 


d'^if/y{x,z,t) 


(12) 
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on  both  material  layers.  Equation  (1 1)  is  the  dilatational  component  and  equation  (12)  is  the 
shear  component  of  the  displacement  field.^  Correspondingly,  the  constants  Cd  and  Cs  are  the 
complex  dilatational  and  shear  wave  speeds,  respectively,  and  are  determined  by 


^  +  2// 
P 


(13) 


and 


(14) 


The  relationships  of  the  Lame  constants  to  the  compressional  and  shear  moduli  are  shown  as 


Z  = 


Ev 

(l  +  i;)(l-2u) 


(15) 


and 


E 

2(1 +  u)’ 


(16) 


2  2 

where  E  is  the  complex  compressional  modulus  (N/m  ),  G  is  the  complex  shear  modulus  (N/m  ), 
and  u  is  the  Poisson’s  ratio  of  the  material  (dimensionless). 


The  conditions  of  infinite  length,  two-dimensional  response  (y  =  0),  and  steady-state 
response  are  now  imposed,  allowing  the  scalar  and  vector  potential  to  be  written  as 

(/>{x,z,t)-(^(z)  exp(iA:^x)  exp(i  cof)  (17) 


and 

y/y{x,z,t)  =  'P(z)exp(iA:^x)exp(i^yO,  (18) 

where  kx  is  the  wavenumber  of  excitation  with  respect  to  the  x-direction  (rad/m),  a  is  the 
frequency  of  excitation  (rad/s),  and  i  is  the  square  root  of-1.  Note  that  equations  (17)  and  (18) 
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are  valid  on  both  layers.  For  acoustic  excitation,  the  wavenumber  is  related  to  the  incident  angle 
by  =  (<y /cy)sin(^) ,  where  6  is  the  incident  angle  (rad)  with  a  value  of  0  corresponding  to 

broadside  excitation. 


Inserting  equation  (17)  into  equation  (11)  yields 


dz^ 


+  a^^(z)^0, 


(19) 


where 


a  = 


(20) 


and 


CO 


Inserting  equation  (18)  into  equation  (12)  yields 


d^'¥(z) 

dz^ 


+/?Mz)  =  0, 


(21) 


(22) 


where 


(23) 


and 


(24) 


The  solution  to  equation  (19)  is 


0(z)  =  A{kj^,o})&x^(\ccz)  +  B{k^,co)s)i.^{-\az), 


(25) 
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and  the  solution  to  equation  (22)  is 


'P(2)  =  C(k^,  <w)exp(i/?z)  +  D(k^,  6;)exp(-i/?z) ,  (26) 

where  B,  C,  and  D  are  constants  that  are  determined  below.  The  displacements  can  now  be 
written  as  functions  of  the  unknown  constants: 

u.{x,z,t)  =  {  \a[A{k^,a)Qxp(iaz)-B{k^,(jo)exxf{-\az)] 

+  iA:^[C(A:^,ry)exp(iyGfe)  +  D(^^,ru)exp(-iy&)]  }exp(i^^x) exp(i<yO  (27) 


and 


{x,  z,t)  =  {  '^k^ {A{k^ ,  co)  exp(i«z)  +  B{k^ ,  co)  exp(-iaz)] 

-\fi{C{k^,co)exp(\fSz)-D{k^,cci)exp{-\Pz)\  } exp(iA:j,x) exp(iru/).  (28) 

The  solutions  of  the  constants  are  determined  by  formulating  the  problem  with  the  presence  of  an 
acoustic  load. 

The  acoustic  pressure  in  the  fluid  is  governed  by  the  wave  equation  and  is  written  in 
Cartesian  coordinates  as"* 


d‘^P^{x,Z,t)  ^^Pa{x,Z,t)  1  d^p^{x,Z,t) 

c^j  dt^ 


(29) 


where  Pg  (x,  z,  t)  is  the  pressure  (  N/m^  ),  z  is  the  spatial  location  (m)  normal  to  the  plate,  and  c j 

is  the  compressional  wavespeed  of  the  fluid  (m/s).  The  acoustic  pressure  is  modeled  as  a 
function  at  definite  wavenumber  and  frequency  as 

p„(x,z,t)  =  Pg(z,k^,G))Qxp(\k^x)ex^(^(ot).  (30) 

Inserting  equation  (30)  into  equation  (29)  and  solving  the  resulting  ordinary  differential  equation 
yields 

Pg(z,k^,(o)  =  G(k^ ,  co)  exp(i;z)  +  Pj  (co)  exp(-i  yz),  (3 1 ) 
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where  the  first  term  on  the  right-hand  side  represents  the  reradiated  pressure  field  and  the  second 
term  represents  the  applied  incident  pressure  field  (the  forcing  function)  acting  on  the  structure. 
In  equation  (31), 


7  = 


(32) 


where  y  is  purely  real  or  imaginary,  depending  on  the  sign  of  the  argument.  When  the  sign  of 
the  argument  is  positive,  the  analysis  is  in  the  acoustic  cone,  and  when  the  sign  of  the  argument 
is  negative,  the  analysis  is  in  the  nonacoustic  region.  For  acoustic  sonar  response,  the  analysis  is 
typically  studied  in  the  acoustic  cone. 

The  transfer  functions  are  now  derived  by  coupling  the  pressure  in  the  fluid  to  the 
displacement  in  the  elastomer.  The  interface  between  the  fluid  and  the  solid  surface  at  z  =  c 
satisfies  the  linear  momentum  equation,  which  is^ 

^  d\.ix,c,t)  _  dp^{x,c,t) 


where  Py  is  the  density  of  the  fluid  (kg/m^).  The  normal  and  tangential  stresses  in  the  system  at 
the  boundary  z  =  c  are 


/  /I  .^.{x,c,t)  .  duJx,c,t) 

r„  {x,  C,  0  =  (/^  +  2//2  )  — -  +  ^2  c,  /) 


dz 


dx 


(34) 


and 


T,,ix,C,t)  = 


du^{x,c,t)  du.{x,c,t) 
dz  dx 


=  0, 


(35) 


where  the  subscript  2  denotes  plate  2,  which  is  the  elastomer.  The  interface  between  the  plates 
requires  four  equations.  The  first  two  equations,  which  are  displaeement  constraints,  are 

u.  ( b,  Olpiate  0|piate  1  (36) 
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and 


(X,  b,  Olpiate  2  Olpiate  1  ’ 


where  plate  1  is  the  metallic  piece.  The  second  two  equations  are  stress  constraints: 


(^»^’0|p,a,e  2  =  (^’^’Olpiate  1 


^■z.(^»^»0|piate2  = 


Finally,  the  normal  and  tangential  stresses  in  the  system  at  the  boundary  z  =  a  are 


^  X  .^,(x,a,t)  .  t9uJx,a,t)  . 

T,,ix,a,t)  =  (^  +  2/2,)  +  ^  =  0 

^  uX 


^  \du^ix,a,t)  ,  du.{x,a,t)'\  ^ 

z.^{x,a,t)  =  /i^  -  =^- 

^  uX 


Combining  equations  (27)  through  (41)  yields  the  eight-by-eight  linear  system  of  equations 


Ax  =  b. 


where,  using  c  =  0,  the  entries  of  equation  (42)  are 


^1,1  ”  “^2^ 


2  ,  Pj^  «2 


A^2=-(X2^-2a2i^2-^K  - 


2  Pj^a^ 


Ai  j  -  —2k^p2Pi  + 


pja?k^ 
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PfCo 

Aa  =  ^ 

y 

(46) 

A.5  - 

(47) 

A,6  - 

(48) 

A,7  - 

(49) 

o' 

II 

(50) 

A,\  =-^M2Ka2^ 

(51) 

A,2  = 

(52) 

^2,3  “  /^2^2  ~  M2^1c^ 

(53) 

A  A  -  f^2P2  ~  fhAl’ 

(54) 

A,5  - 

(55) 

A, 6  ~ 

(56) 

Aj  ~ 

(57) 

A, 8  ~ 

(58) 

^3,1  =0«2)exp(ia2^), 

(59) 

^3,2  =(-i«2)exp(-ia2^)), 

(60) 

A,3  =  (iAh^pOA^)’ 

(61) 
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4.4  =(i^;.)eXP(-jA^X  (62) 

4.5  =  (-  *<^i  )exp(ia]^>),  (63) 

4.6  =  0^1  )exp(-iai6),  (64) 

4.7  =  (-  i4)exp(iA^>),  (65) 

4.8  =  (-  i4  )exp(-iA^)»  (66) 

4.1  =04)®^P(*^2*X  (67) 

4.2  =04)®xp(“i^2^)’  (6^) 

4.3  =  (-  iy^2  )exp(i/?26),  (69) 

4.4  =  0  A  )exp(-iy^2^)’  (76) 

^^4  5  =  (-  \k^  )exp(i  a^b),  (7 1 ) 

A^f,  =  (-  iA:^)exp(-ia,6),  (72) 

4.7  =(iA)exp(iAi*),  (73)- 

4.8  =  (-  iA  )exp(-iA6),  (74) 

4.1  =  (~ <^24  “ 2a2>“2  ~ 4^1  )exp(i«2^)»  (75) 

4.2  =  (“<^24  “2a2/^2  ~4^^)®^P(“j^2^)’  (76) 

4.3  =  (-  '^KPiMi  )exp(i^2^)’  (77) 

4.4  =  )exp(-iA6),  (78) 
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^5,5  =  +  2«]Vi  +  )e\p(io'jb),  (79) 

^5,6  =  +  2«iVi  +  )exp(-iajb),  (80) 

^5,7  =  (81) 

As  =  (-  )exp(-i/?]6),  (82) 

A.1  =  (-  2/y2A:^«2)exp(ia26),  (83) 

A,2  =  )exp(-i«26),  (84) 

A, 3  =  [^^2^2  -  M2^1  )exp(i/?26),  (85) 

Aa  =  (^2/^2  -  /^2^^  )exp(-i;^2^)’  (86) 

Aa  =  (2  A  )exp(i  a^b),  (87) 

A,6  =  (-  2a  )exp(-iori6),  (88) 

^6,7  =  (-  M\PI  +  )exp(iA^),  (89) 

A,%  =  (-  A  +  A^x  )exp(-iA^),  (90) 

A,i=0,  (91) 

A,2=0,  (92) 

A,3=0,  (93) 

A,4=0,  (94) 
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Ayj  =  (-  ajAj  -  2aj  ]exp(iQ',o), 

(95) 

A,6  “  (“  A  “  )exp(-iaia). 

(96) 

^7  7  =  (-  )Qxp(SP^d), 

(97) 

^7,8  =  (2A:;,A/“i  )exp(-i^,fl). 

(98) 

o' 

11 

(99) 

> 

po 

K> 

II 

o 

(100) 

r\ 

o 

II 

m 

(101) 

> 

00 

4S» 

II 

O 

(102) 

^8  5  =  (-  2M\Ka\  )exp(ia,fl), 

(103) 

^8,6  =(2/^i^;(«i)exp(-iQ;,a), 

(104) 

^7  =  U  )exp(iA«). 

(105) 

4,8  =  (a A?  -M\^l  )exp(-i/?,a). 

(106) 

Xi,  =A2ik^,co), 

(107) 

Xii  =  B2{kj^,cd), 

(108) 

X3]  =C2(A:^,<w), 

(109) 

X4,  =Z)2(A^,<y), 

(110) 

^5,1  = 

(111) 
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(112) 

xi,x=c^(.KM,  (113) 

H,x=Dxik^,o)),  (114) 

b,,=-2P,{o)\  (115) 

62,1-0,  (116) 

63,1=0,  (117) 

64,1=0,  (118) 

65,1=0,  (119) 

65,1=0,  (120) 

67,1=0,  (121) 

and 

68,1=0.  (122) 


Using  equations  (43)  through  (122),  the  solution  to  the  constants  C, ,  £>j ,  A2,  B2,  C2 , 

and  D2  can  be  found  by 

x  =  A"’b.  (123) 

The  transfer  function  between  the  normal  displacement  at  location  z,  when  z<b,  and  the 
incident  pressure  (acoustic  load)  can  be  written  as 
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"  Pi(oy) 


=  (k^ ,o))iai  exp(ia,r)  -  (k^ ,  cd)  exp(-iaiz) 

+  C^{k^,co)\k^  exp(i>0j z)  +  D^{k^, od) ik^  exp(-iy0i z)  (124) 


or,  when  z  >  6 ,  as 


T{z,k^,od) 


^  U.(z,k^,co) 

P,ico) 

=  A2{k^,od)\a2  exp(i «2 ^) “ ^2 (A:^ , 6?) i 0^2  exp(-i ^2 

+  C2 {k^ ,  cd)  \k^  exp(i>92^)  +  P>2 i^x  >  exp(-i/?22> 


(125) 


The  transfer  function  between  the  tangential  displacement  at  location  z,  when  z<b,  and  the 
incident  pressure  can  be  written  as 


T(z,k^,o)) 


^Uxiz,k^,o)) 

PjiCO) 

=  A^  (k^ ,  cd)  ik^  exp(ia,z)  +  Bi  (k^ ,  cd)  i^^  exp(-i a^z) 

+  -C,  (k^ ,  cd)  \jd^  exp(iy5iz)  +  (k^ ,  co)  i/9,  exp(-iy9iz) 


or,  when  z  >  6 ,  as 


T(z,k^,6)) 


Ux{z,kx,co) 

Pi{k^,co) 

-  A2  (k^ ,  cd)  ik^  exp(ia2z)  +  B2  (k^ ,  cd)  \k^  exp(-i otj^) 

+  -C2 {k^ ,  od)  i/92  exp(iy92z)  +  £>2 {k^ ,  od)  1P2  exp(-i^2^)- 


(126) 


(127) 
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The  transfer  function  between  the  normal  stress  at  location  z,  when  z<b,  and  the  incident 
pressure  can  be  written  as 


T{z,K,o)) 


S„{z,k^,(o) 
Pi  (CO) 


=  A^  (k^ ,  <w)(-  -  la^jUi  -\kl  )exp(i a^z) 

B^{k^,co^  a\  \-2a\  -  lykl^&xip{-\a^z) 

+  C,  (k^ ,  oj\-lkJ^iu^ )  exp(iy5,z)  +  (k^ ,  co\2k^p^iu^ )  exp(-/y0,z) 


(128) 


or,  when  z  >  6 ,  as 


T(z,k^,co) 


_  S..(z,k^,(o) 

P,(co) 

=  A2  (k^ ,  <y)(-  «2  ^  -  2^2  )exp(i  ^2^) 

+  B2(k^,Q)i^ ^2 - 2al/J2  - ^kl )exp(-iQr2z) 

+  C2  ,  (o\-2kj2Bi )  exp(iy?2z)  +  A  (^x  >  <yX2^j,  A>"2 )  exp(-i 


(129) 


The  transfer  function  between  the  shear  stress  at  location  z,  when  z<h,  and  the  incident 
pressure  can  be  written  as 


T(z,k^,o)) 


S^(z,k,,a)) 

PliK’CO) 


=  -Ai  (kj^ ,  co)  (2//i^^cir,  )exp(i  a^z)  +  5,  (k^ ,  o)){2/u^k^ai  )exp(-i  qtjz) 

+  Ci(k„co)[^i(^i  -A^)]exp(i/?iz)  + A(^;c»^)[m(A^  -^^)]exp(-iAz) 


(130) 


or,  when  z>b,  as 


T{z,k^,cd) 


S._^{z,k^,(o) 

PjiK^co) 


=  -^2(^;c’^X2/^2^;r«2)exp(ia2^)  +  A(^;c>^)(2/^2^;t«2)eXp(-iQr2Z) 

+  C2(k^,0})[^2(P2  -^x)]exp(iA^)  +  A(^x>^)[/^2(>^2  - ^^)]exp(-i/?24 


(131) 
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Equations  (124)  through  (131)  are  the  open-form  solutions  to  the  problem,  and  these  equations 
provide  the  correct  answer.  However,  the  physics  are  not  evident,  and  solving  a  corresponding 
inverse  problem  is  typically  problematic. 
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3.  REDUCED-ORDER  MODEL 


A  reduced-order  model  presents  the  transfer  functions  in  a  closed  form,  which  is  a  more 
desirable  expression  than  the  open-form  solutions  derived  in  section  2.  Deriving  a  reduced-order 
model  is  facilitated  by  changing  the  model  of  the  bottom  (metallic)  plate  from  an  elastic  medium 
to  a  continuous  mass  that  is  infinitely  stiff.  This  is  a  reasonable  assumption  because  at  high 
frequencies  the  metallic  plate  appears  rigid  to  acoustic  plane  waves,  and  flexural  wave  energy 
does  not  have  a  significant  dynamic  effect  on  the  plate  response.  Although  this  assumption 
limits  the  model,  it  is  reasonably  accurate  in  the  acoustic  domain.  The  boundary  conditions  of 
the  model  at  z  =  c  remain  equations  (33),  (34),  and  (35).  The  boundary  conditions  at  z  =  6 
become 


(132) 


and 


T,^{x,b,t)  = 


^^{x,b,t)  ^.(x,b,t) 
dz  dx 


=  -M 


d^uj^x,b,t) 


(133) 


where  M  is  the  mass  per  unit  area  (  kg/m^ )  of  the  metallic  plate.  For  brevity,  the  subscript  2  is 
now  dropped  from  the  analysis  and  unsubscripted  material  constants  are  related  to  the 
elastomeric  plate.  Once  these  boundary  conditions  are  applied  to  equations  (27)  and  (28),  the 
result  is  a  four-by-four  linear  system  of  equations  that  can  be  written  as 

Ax  =  b,  (134) 

where  the  entries  of  equation  (134)  are 

2 

^  p  fCO  cx 

A^^=-a^X-2a^/u-Xkl-\-^ - ,  (135) 

y 

A^-^=-a^X-la^  (136) 

r 
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4.3=  2M/^+ 

(137) 

44=2M/^+  ^ 

(138) 

7^2  ,  = 

(139) 

yi2  2  =  2fikj^a, 

(140) 

4.3  =  mP^-m^I^ 

(141) 

^2,4  = 

(142) 

^3 1  -  (“  “  2a^/2  -  /lA:^  +  i^y^«M  )exp(iorZ)), 

(143) 

4.1  =  - 2a^ii - Xkl  -loPaM^Qxpi-xab), 

(144) 

^3  3  “  (“  ^^xPm  +  i<y^A:^M)exp(i/?&), 

(145) 

^3,4  =  ^^xPm  +  'xaPk^M^ex.p(-\pb), 

(146) 

=(-2//A^Gr-i«^A^M)exp(i«Z>), 

(147) 

A^  2  -  i^pkx^  -  ’^oPk^M^QXTp{-iab), 

(148) 

4.3  =  [mP^  -  Juki +\oPpM)Qx^(\Pb), 

(149) 

A4  4  =  {fiP^  -M^l  -\oP pM)Qxp(-\Pb), 

(150) 

\\  =  4K^o)), 

(151) 

^2,1  = 

(152) 
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Xu^C{k^,o}), 


(153) 


=  D(k^,o)), 

(154) 

(N 

1 

II 

(155) 

0 

II 

0 

(156) 

II 

0 

(157) 

and 

641=0.  (158) 

Using  equations  (135)  through  (158),  the  solution  to  the  constants  A,  B ,  C ,  and  D  can  be 
found  by 

x  =  A"*b.  (159) 

The  right-hand  side  of  equation  (159)  can  be  solved  in  a  closed-form  solution  by  inverting 
the  A  matrix  and  multiplying  the  results  by  the  b  vector.  The  constants  in  the  x  vector  are  known 
explicitly  and  can  be  re-inserted  into  equations  (125),  (127),  (129),  and  (131)  to  solve  for  the 
displacements  and  stresses.  In  addition,  the  determinant  of  the  system  matrix  is  available  as 

det(A)  =  4p^A{k^,G)),  (160) 

where 

A{k^,co)  =  di+d2  cos(ah) cos{ph)  +  d^  cos{ah) s\n{fih) 

+  d^  sm{ah)  cos(y56)  +  d^  sin(ah)  sin(yS6),  (161) 


with 


d,  =4apk^^(^^  -k^j[2(P^  -k^J-{M / p)\/3^  +k^J^], 


(162) 
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=  -a/3[ul{P^-klfHMIp)\P^+kl)\5kt-2klp^+P^) 

+  \{pj  lpy){Mlp){/3^+klf\ 

d,  =  m{/3^+klf[ipj !  pyMP^ -ip^l^kl) 

+  -\{Mlp){p^-6p^kl  +  kl)  +  {pj  lpy){Mlpfkl{P^+klf\  (164) 

=  xpip'^  +klf[4{pj  I  py)a^kl 
+  i(M  /  p\p'^  +  2aJt^  -  kl  ){p'^  -  2ak^  -  k^ ) 

+  ipy /py)(M/pfa\p^  (165) 


and 

rf5=(A"-*?)‘'+16*>V  + 

(M  /  +(,/}^ -klf] 

+  -\(pjl  py\MI  p)a^(f  ^kl)‘ .  (1««) 

Once  the  determinant  of  the  system  is  known,  the  transfer  function  between  the  tangential 
displacement  and  the  incident  pressure  can  be  computed  by 

nz,k^,(o)  =  =  Ul{z,k,,(o) 

'  ’  '  Pjico)  p^ik^,co) 

where 


U I (z, k^,o))  =  ik^Ci  cos(az)  +  iA:^C2  sin[a(/j  +  z)] sin{ph) 

+  iA:^C3  cos[«(/z  +  z)]  sin(Ph)  +  ik^c^  sin[a(/2  +  z)]  cos(ph) 

+  xk^c^  cos[a{h  +  z)]cos(y9/z)  +  iy&r^  cos(y9z) 

+  xPc-j  s,xrx[P{h  +  z)]sin(af^)  +  iyficg  cos[yff(/i  +  z)]sin(cir/?) 

+  xPcg  sin[y0(/?  +  z)]cos(Qr^)  +  iy^Cjo  cos[y9(^  +  z)\cos{ah),  (168) 
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with 


c,  =‘iapkl[^0^ -kl)-(Ulp)\p^*klf\ 

(169) 

-kljiP^ -klf  +(,Mlpykl(p^  +*,")"], 

(170) 

cj  =  2a(MlpXP" -klXP^+klf, 

(171) 

c,  =-2piMlptp'^-kl\P^^kl)\ 

(172) 

c,  =-2ap(JI^ -kl'^kl*(M  !  p)\p^  +^^1 

(173) 

c,=-2cdi,(p^ -k]i2(P^ -kl)-^M  !  p)\p^  *klf\ 

(174) 

Ct  =-ia^pkjfikl+{M !  pf(J3Uklf\ 

(175) 

Ct=-*aH,(Mlp)kP^+kl)\ 

(176) 

c.^'kapk.iMIpXP^^klf, 

(177) 

and 

<'10  -klf  +(M  !  pfkl(P^  +klf\  (178) 

The  transfer  function  between  the  normal  displacement  and  the  incident  pressure  can  be 
computed  by 


T(z,k^,0)) 


U,(z,k,,o))  _Uliz,k,,co) 
P,(o]i)  ^{k^,o)) 


(179) 


where 
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U I (z,k^,o))  =  -ac,  sin(Qs)  +  ac2  cos[Qr(/z  +  z)] sm(/3h) 

+  -ac^  sin[a(//  +  z)]sin(y6/z)  +  ac^  cos[ar(/z  +  z)]  cos(/?/?) 

+  -ac^  sin[a(/7  +  z)]cos(j3h)  +  sin(y&) 

+  -k^C’j  cos[j3(h  +  z)] sin(a/z)  +  k^c^  sin[y0(/z  +  z)] sin(a/?) 

+  -k^CgCos[fi(h  +  zy\cos(ah)  +  k^c^Qsin[P(h  +  z)]cos(ah).  (180) 

The  transfer  function  between  the  shear  stress  and  the  incident  pressure  can  be  calculated  using 


T(z,k„o)) 


^zx  (^’  _  ^:x  (•^  ? 

Pj(co)  Mk^,o)) 


(181) 


where 


{z,k^,co)  =  -2i ak^Cy  sin(az)  +  2i ak^C2  cos[a{h  +  z)] sin{ph) 

+  -liak^c^  sin[«(/j  +  z)]sin(y9/?)  +  2\ak^c^  cos[a(h  +  z)]cos{/Sh) 

+  -Ixak^c^  sin[«(^  +  z)]cos(y6/7)  -  i(/?^  -  kl)C(^  sin(y0z) 

+  i{P^  -  kl )c-j  cos[y9(/?  +  z)] sm{ah)  -  \{P^  -kl)c^  sin[yff(/z  +  z)] sinfor/z) 

+  \{P^  -kl)cgCos[p{h  +  z)]cos{ah)-\{P^  -kl)CyQS\n[P{h  +  z)']cos{ah).  (182) 


The  transfer  function  between  the  normal  stress  and  the  incident  pressure  can  be  calculated  using 


S,,(z,k^,co)  _  Sl,{z,k^,co) 
Pjico)  A(k^,co) 


(183) 


where 

Sl.{z,k^,co)  =  -{P^  -kl)cy  cos(az)-(/?^  -^1)^2  sin[a(/z  +  z)]sin(/?//) 

+  -{P^ -^1)^3  cos{a{h  +  z)]s.m{pK)  +  -{p^ -kl)c^  sin[«(/z  +  z)]cos(/%) 

+  -{P^  -  k] )c^  cos[ar(/;  +  z)]  co^{ph)  +  ipk^pf^  cos(>9z) 

+  2pk^C2  sin[/9(/?  +  z)]sin(«/z)  +  2pkj^c^  cos[P{h  +  z)]sin(cr//) 

+  2pk^c^  sm[P{h  +  z)]cos(a/2)  +  2pk^Cy^  cos[P(h  +  z)]co^{ah).  (184) 
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The  reflected  pressure  field  in  the  fluid  is  equal  to 


P,t(zy,k^,co)  = 


ir 


U._{b,k^,(o)  +  \ 


exp(ifZy), 


(185) 


with  Zj  being  the  position  where  the  field  is  evaluated  (m).  The  total  pressure  field  is  the  sum  of 
the  reflected  pressure  field  and  the  spatially  phase-shifted  excitation  level.  It  is  written  as 


Pfoiai  (^/>k^,0})  =  Pjfizy,k^,0))  +  Pj  (0))  exp(-iyzy). 


(186) 


Equations  (132)  to  (186)  are  the  reduced-order  model. 
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4.  COMPARISON  OF  MODELS 


By  comparing  the  full-  and  the  reduced-order  models,  the  effects  of  modeling  the  metallic 
backing  piece  as  a  discrete  mass  rather  than  as  a  thick  plate  can  be  discerned.  The  baseline 

model  includes  an  elastomer  with  the  properties  of  /?  =  1200 kg/m^,  X  =  9.8e8(l-0.01630N/m^, 

ju  =  le7(l  -  0.2i)  N/m^ ,  and  h  =  0.0508  m;  a  steel  backing  plate  with  the  properties  of 

p  =  7850  kg/m^ ,  ;i  =  1.21el  1(1-  0.01i)N/m^ ,  p  =  8.04el  0(1  -0.01/)  N/m^ ,  and  //  =  0.0127  m; 

and  fluid  with  the  properties  of  p  =  1 025  kg/m^  and  =  1 500  m/s .  The  model  is  evaluated  at 

z  =  -0.0381  m.  Figures  2  through  1 1  present  the  system  responses  versus  incident  angles  at 
frequencies  from  10  kHz  to  100  kHz  in  10-kHz  increments.  In  each  of  the  figures,  the  solid  line 
is  the  double-plate  model  developed  in  section  2  and  the  “x”  marks  are  the  reduced-order  model 
developed  in  section  3.  In  addition,  the  two  plots  in  the  upper  left  of  each  of  the  figures  are  the 
magnitude  and  phase  angle  of  the  normal  displacement  divided  by  incident  angle  pressure  versus 
incident  angle.  The  two  plots  in  the  upper  right  are  the  magnitude  and  phase  angle  of  the 
tangential  displacement  divided  by  incident  angle  pressure  versus  incident  angle.  The  two  plots 
in  the  lower  left  are  the  magnitude  and  phase  angle  of  the  normal  stress  divided  by  incident  angle 
pressure  versus  incident  angle.  The  two  plots  in  the  lower  right  are  the  magnitude  and  phase 
angle  of  the  shear  stress  divided  by  incident  angle  pressure  versus  incident  angle.  It  is  evident 
that  at  most  frequencies  the  reduced-order  model  is  equivalent  to  the  double-plate  model.  The 
two  models  diverge  somewhat  at  some  specific  incident  angles.  This  disagreement  corresponds 
to  elastic  waves  in  the  full-order  model  that  are  not  present  in  the  reduced-order  model.  If  the 
mass  of  the  metallic  plate  is  increased,  the  comparison  will  be  more  favorable.  If  the  mass  of  the 
metallic  plate  is  decreased,  the  two  models  will  begin  to  deviate  from  one  another. 
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Figure  8.  System  Response  Versus  Incident  Angle  for /=  70,000  Hz 
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5.  CONCLUSIONS 


This  report  has  developed  a  reduced-order  model  that  provides  a  closed-form  solution  to  the 
problem  of  a  fluid-loaded,  elastomeric-coated  metallic  plate.  This  solution  form  is  useful 
because  the  dynamics  of  the  problem  are  now  explicit  and  the  contribution  of  each  term  to  the 
solution  is  clearly  evident.  This  solution  has  been  compared  to  the  open-form  (implicit)  solution 
of  a  fluid-loaded,  double-thick  plate  system  from  10  kHz  to  100  kHz.  The  result  is  that  the 
reduced-order  model  is  almost  identical  to  the  open-form  solution  for  all  the  field  variables.  This 
solution  is  limited  to  the  case  where  a  stiff  plate  is  covered  by  a  relatively  soft  plate.  Further 
work  in  this  area  should  be  undertaken  with  experiments  so  that  actual  data  can  be  compared  to 
the  models. 
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